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r-j • Abstract. We discuss the propagation of spectral line and continuum radiation in a 

clumpy medium and give general expressions for the observed absorption or emission 
*vj ■ from a cloud population. We show that the affect of the medium dumpiness can 

^\1 ' usually be characterised by a single number multiplying the mean column opacity. 

Our result provides a simpler proof and generalization of the result of Martin et al. 
(1984). The formalism provides a simple way to understand the effects of clumping 
^ ' on molecular line profiles and ratios; for example how clumping effects the interpre- 

^^ , tation of ^^CO(l-O) to ^■^CO(l-O) line ratios. It also can be used as a propagation 

Qy ' operator in physical models of clumpy media where the incident radiation effects the 

spectral line emissivity. We are working to extend the formalism to the propagation 
of masers in a clumpy medium, but in this case there are special difficulties because 
formal expectation values are not characteristic of observations because they are 



^— V , biased by rare events. 

r^^ • Keywords: molecular lines, spectra 

o 

52 ' 1. Introduction 



The problem of the propagation of radiation in a non-uniform medium 
is a general one in astrophysics, which arises in many different contexts. 



C^ ■ We consider here the common case that the medium consists of clumps 

or clouds with a low volume filling factor, but potentially high covering 
factor. These clumps or clouds may possess a range of properties. This 
paper gives a summary of our derivation and results, a more detailed 
proof will be published elsewhere (Conway, Elitzur and Parra 2004, in 
prep). 

The specific problem we address is the expected response of a tele- 
scope beam observing a low filling factor clumpy region where the 
beamwidth is much larger than the typical cloud diameter. In this case 
we can show that that expectation value of the beam average equals 
the ensemble expectation value over many realisations of observations 
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along the line of sight (LOS). We first (see Section 2) analyse the case 
of continuum (e.g. dust) clouds absorbing background radiation, then 
spectral line absorption (Section 3). In Section 4 we consider spectral 
line emission and apply the formalism to the interpretation of CO line 
ratios. Finally in Section 5 we draw conclusions and describe future 
work. 



2. Continuum Absorption Radiation Transfer 

Let us consider continuum radiation transfer in a clumpy medium of 
dust clouds. First assume that all clouds/clumps are identical. Assume 
further that radiation / entering a clump emerges as I ■ R, where the 
clump response is R. In terms of cloud opacity R = e~'^. Hitting suc- 
cessively k clumps that all have the same response factor produces the 
aggregate response R^. Implicit in this assumption is that the interac- 
tion with the radiation does not change the clump response factor R. 
It can be shown that provided the clump volume filling factor is small 
the probability of getting k clumps along a given line of sight always 
obeys a Poisson distribution, characterised by one number, the mean 
N along a LOS. In this case of one clump type it can be shown that 
the mean response is 

R = e-\ t = N{l-e-^) (1) 

as previously derived by Natta & Panagia (1992). 

Next consider the case of many types of clumps with response factors 
Ri (these could be clump types having different optical depths) each 
obeying Poisson statistics. If the fraction of clumps of each type is fi 
then (see Conway, Elitzur and Parra 2004, in prep) 

i? = exp (iV[i?av - 1]) , where R^, = Y,hRi (2) 

i 

This means that the problem of a mixture of cloud optical depths (tj) 
is equivalent to that of the single cloud type having an average Tav 
obtained from 

i?av = e-"-=^/,e-"^ (3) 

i 

We can generalize the result in Eq. (2) to the case where the clump 
type varies continuously. For instance if the clump type is described 
by continuous variables t and r (which may correspond for instance to 
peak optical depth and distance from the centre of clouds which have 
spatial structure) then we have 



conway.tex; 2/02/2008; 11:42; p. 2 



Spectral lines in a clumpy medium 



R^^= / R{r,t)7]{r,t)drdt = / e-^^'''^^r]{r,t)drdt (4) 



and likewise for any number of dimensions of continuous parameters, 
where r] is the distribution function of clump types. 



3. Line Radiation Transfer - Absorption 

In the spectral line absorption problem we aim to calculate the ex- 
pectation value of a frequency dependent absorption response function 
given a population of clouds. This cloud population has in general a 
distribution of central LOS velocities and a range of spectral profiles 
and spatial structures. Observed absorption lines are generally pub- 
lished as a function of Doppler velocity v, the response at velocity 
R{v) being understood to represent the response at observing frequency 
z^ = 1^0(1 + v/c) where z^o is the line's rest frequency. Let us concentrate 
on calculating the response at a particular velocity v, R{v). In principle 
clouds having all possible centre velocities v' contribute to the spectral 
response at velocity v. The centre velocity of each cloud v' can simply 
be thought of as yet another continuous variable describing the cloud 
properties, so the full distribution function becomes r]{r, t, v'). For spec- 
tral line clouds the response function is Rv{r, t, v') which if cloud types 
are defined to have fixed spectral profiles which don't vary with centre 
velocity can be written R{r,t,v — v'). Hence applying Eq. (2) and (4) 
we obtain 

R{v) = exp{N[R,,{v)-l]) (5) 



where Rav{v) = / / R{r,t,v — v')'r]{r,t,v')dv'drdt (6) 

As before N is the mean number of clouds (at any velocity) per LOS. 
Raviv) can be interpreted at each r,t value as a convolution of the 
cloud velocity profile for this r,t, R{r, t, v) with the velocity distribution 
function of the cloud population rj{r, t, v) followed by an integration 
over all r,t (and any other continuous cloud parameters). 

If the fractional distribution of clouds of different kinds is indepen- 
dent of velocity, then 7]{r,t,v) = f{r,t)q(v) is a separable function. In 
this case we can define a mean cloud spectral response 

Rs{v)= f f R{r,t,v)f{r,t)drdt (7) 
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In this separable case the expectation value of the emerging spectrum 
is 

R{v)=exj>{-X{v)*C{v)) (8) 

In this expression * denotes convolution, A(f ) = Nq{v) describes the 
cloud velocity distribution; more precisely defined as the mean number 
of clouds per unit velocity per LOS. The quantity C{v) = 1 — Rs{v) = 
< l — e~'^^'"> >, where t{v) is the opacity profile of each clump type, and 
the brackets denote an average over the clump population. The C{v) 
can be interpreted as the mean cloud emission profile. In the normal 
case that the X{v) distribution is much broader in velocity than C{v) 
we obtain the good approximation that R{v) ~ exp (— CyiA(f )) where 
Ca = / C{v)dv has dimensions of velocity and can be considered the 
equivalent width of the mean cloud profile. 

Another way to view the effect of clumping is to note that 

Riv) ^ exp i-Kftotiv)) (9) 

where ftot{v) = X{v)* < t{v) > is the mean line of sight opacity and 
i^ is a factor taking into account clumping which is 

J < T(V) > dv 

Remarkably we find that the effect of an arbitrarily complex clump 
distribution can usually be reduced to a single number, K. This factor 
K reduces the effective opacity of a clumpy medium compared to a 
smooth gas of the same mean column density. Martin et al. (1984) 
partially derived our results, but their derivation was less direct and 
the full generality of the result was not stated. 



4. Emission Line Spectra 

The expressions for opacity derived above can be adapted (Conway, 
Elitzur and Parra 2004, in prep) to predict the emissivity in the case of 
a population of clouds with a mixture of excitation temperatures and 
statistical properties along the line of sight. Here we discuss only the 
simplest case that all clouds have the same excitation temperature Tg. 
In this case the observed spectral brightness temperature is 

T{v) = T,[l - R{v)] =T,[l- exp {-Kftot{v))] (11) 

To illustrate the impact of clumping on observed line profiles and line 
ratios we show an example in Fig 1 involving the ^^C0(l-O)and ^^C0(1- 
0) lines. In this example the cloud population consists of only one clump 
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type with a Gaussian opacity velocity profile t{v), and a much broader 
Gaussian X(v) cloud velocity distribution. The expectation value of the 
total opacity along each LOS ftot{v) = A(w) * t{v) is therefore also a 
Gaussian whose peak we term the 'peak total opacity'. 

The different panels for Fig 1 have different combinations of the 
^^CO(l-O) peak total and peak cloud opacities. The ratio of intrinsic 
^^CO(l-O) and ^^CO(l-O) opacities is set to 60 so that in all the panels 
of Fig 1, the ^^ CO (1-0) total and cloud opacities are << 1 and hence 
are optically thin. In all cases therefore the ^^ CO (1-0) line profiles are 
Gaussian with peak value (in units of Tg) equal to the total ^^CO(l-O) 
opacity divided by 60. In contrast the ^^ CO (1-0) line profiles and their 
ratio to the ^^CO(l-O) lines depend critically on the clump opacity 
which affects K. Only for the case in the top left panel, which has both 
a low total and low clump opacity, is the observed line ratio equal to 
the intrinsic value. If either or both total and clump opacities are larger 
than 1 then a smaller line ratio is observed, but it is impossible just 
from the line ratio at line centre to distinguish the different cases. 

An important consequence of Eq. (9) is that when the cloud velocity 
distribution is much broader than individual clouds it is in princi- 
ple impossible from the ^^CO(l-O) and ^^CO(l-O) spectra alone to 
separate the effects of clumping and variations in the intrinsic ^^C0(1- 
0)/^^CO(1-0) opacity ratio. In going from panels on the top row with 
low opacity clumps to the bottom row with large clumps the ^^C0(1- 
0) clump factor K is reduced. The resulting spectra are the same as if 
there was no clumping but the intrinsic ^^CO(1-0)/^^CO(1-0) opacity 
ratio was reduced. 



5. Conclusions and Future Work 

We have derived general expressions for both continuum and spectral 
line radiation propagation in a clumpy medium. These expressions hold 
provided only that the clumps have a small volume filling factor and 
that the local emissivity is not itself affected by the local line spectra. 
Eq. (8), (9) and (10) provide simple expressions that can be used by ob- 
servers to interpret observed line profiles. In Section 4 we illustrated the 
application of this equation to CO line emission profiles from a clumpy 
medium. Although the results presented in Section 4 are not new our 
equations provide a way to understand line profiles from clumpy media 
in a systematic and quantitative way. 

We are working to extend our formalism to predict variances and 
higher order statistics of spectra from clumpy media. We also hope to 
use our formalism as a propagator term in physical spectral line models 
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Figure 1. Example ^■^CO(l-O) and ^^CO(l-O) line profiles (shown respectively by 
solid and dashed lines) for different combinations of cloud and mean column opaci- 
ties. To make them more visible the ^■'CO(l-O) profiles have been multiplied by 10. 
The intrinsic ratio of opacity in the two transitions is assumed to be 60. Both the 
internal cloud opacity profiles of the clouds in the two transitions and the cloud 
velocity distribution are assumed to be Gaussian with the latter having a FWHM 
10 times larger. The top row shows the case of low opacity in each ^^CO(l-O) cloud, 
and a varying mean total ^^ CO (1-0) opacity. The bottom row shows the case of 
high opacity ^^CO(l-O) clouds. Note the different vertical scales of each plot. Only 
when the total and cloud ^^ CO (1-0) opacities are both less than one is the observed 
^^CO(l-O) to ^^CO(l-O) ratio close to the intrinsic value of 60. 

of clumpy media. In general in such models emissivity will depend on 
the local spectrum so such models must be iterative. Finally we are 
actively working (see Parra et al. these proceedings) on the problem 
of maser propagation in a clumpy medium. Formally exactly the same 
equations apply as in the thermal case, the only difference is the use of 
negative opacities. However as well as physical complications of maser 
saturation there is a more subtle effect that even for unsaturated masers 
the formal means predicted by the equations are not good predictors 
of observations. The reason for this is that these mean values are dom- 
inated by events causing very large amplification which are also very 
improbable, such that they never occur within a typical source. 
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